New Results From Lattice QCD: Non-Perturbative Renormalization and 
Quark Masses. 

V. Gimenez a * 

a Dep. de Fisica Teorica and IFIC, Univ. de Valencia, 
Dr. Moliner 50, E-46100, Burjassot, Valencia, Spain. 

For the first time, we compute non-perturbatively, i.e. without lattice perturbation theory, the renormalization 
constants of two-fermion operators in the quenched approximation at /3 = 6.0, 6.2 and 6.4 using the Wilson and 
the tree-level improved SW-Clover actions. We apply these renormalization constants to fully non-perturbatively 
estimate quark masses in the MS scheme from lattice simulations of both the hadron spectrum and the Axial 
Ward Identity in the quenched approximation. Some very preliminary unquenched Wilson results obtained from 
the gluon configurations generated by the T^L Collaboration at j3 = 5.6 and Nf = 2 are also discussed. 



1. Introduction. 

The values of the quark masses are very impor- 
tant for phenomenology. Quark masses are free 
parameters of the QCD Lagrangian which can- 
not be determined within QCD from theoretical 
considerations only. Since quarks are confined in- 
side hadrons, they are not observed as physical 
particles. Therefore, quark masses cannot be 
measured directly and their values depend 
on how they are defined from observations, 
i.e. they are scheme and scale dependent 
quantities. 

The QCD Lagragian has a chiral symmetry for 
vanishing quark masses that is spontaneously bro- 
ken by dynamics and explicitly broken by quark 
masses. The masses of the Goldstone bosons it, 
K and r\ can be computed in a systematic ex- 
pansion in the quark masses, the so-called chi- 
ral perturbation theory (%PT), plus certain un- 
known non-perturbative parameters, A, B, C, so 
that to lowest order one has: 

m PS — C (mi +TO2) 
my = A + B(m 1 +m 2 ) (1) 

where mps (my) is the mass of a pseudoscalar 
(vector) meson in terms of the quark masses m\^- 
As a consequence, xPT can only estimate ratios 
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of quarks masses but not their absolute val- 
ues. Therefore, to determine the absolute nor- 
malization of quark masses, we have to use meth- 
ods which go beyond X-PT, such as QCD Sum 
Rules or Lattice QCD. In this work, we will use 
Lattice QCD to compute the meson masses and 
then, by utilizing the xPT results, extract the 
quark masses. As an alternative procedure, we 
will also simulate the Axial Ward Identity on the 
lattice to calculate quark masses. We do not com- 
pute the bottom quark mass which can be ex- 
tracted from the Lattice HQET Q . 

2. The lattice definition of quark masses. 

In this work, we will consider the Wilson and 
the tree-level improved SW-Clover lattice actions 
in Euclidean space. They are parameterized in 
terms of (/3, re) which are the only tunable input 
parameters related to the bare lattice coupling 
constant, go, and the bare lattice quark mass, m , 
via (3 = 6/ <7q and re = 1/(2 m^a + 8) respectively, 
with a the lattice spacing. As is well known, these 
actions contain a term proportional to a, the so- 
called Wilson term, which is needed to avoid the 
fermion doubling problem, but explictly breaks 
the SU(Nf) L x SU(N f ) R chiral symmetry. On 
the lattice, the chiral properties of QCD are lost 
so that the chiral limit does not correspond to 
vanishing bare lattice quark masses and operators 
of different chirality can mix among themselves. 
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In spite of this, it is possible to find some lattice 
currents which are partially conserved and obey 
the Current Algebra in the limit a — > || . 
Vector Ward Identities (VWI) are obtained 
by applying infinitesimal vector flavour transfor- 
mations to the quark fields of the lattice actions 
and exploting the invariance of the path integral 
under local changes of the fermionic variables. 
Using them, it is easy to show that there is a lat- 
tice vector current, V^, which is both conserved 
in the limit of degenerate quarks, i.e. Z v — 1, and 
improved, i.e. its matrix elements have errors of 
order 0(<7ga) only. For the Wilson action, the lat- 
tice VWI between on-shell hadronic states a and 
f3 takes the form 
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where K% t 2 are the hopping parameter of the 
quarks of the hadron, S = g~2<7i is the bare scalar 
density and is the forward lattice derivative. 
By comparing with the continuum QCD VWI, 
one obtains the relation between the renormalized 
lattice quark mass and the hopping parameter 



mpa = A s ma = Z s —I 

2 \ Hi Kif- 



(3) 



where n c is the critical value of the hopping pa- 
rameter corresponding to a vanishing quark mass, 
and Z$ is the renormalization constant of the 
scalar density. For the SW-Clover action, a sim- 
ilar relation can be derived (see ||). In order 
to determine k, one could compute the ratio of 
the vector and scalar densities on the lattice but, 
unfortunately, the scalar density matrix elements 
turn out to be extremely noisy, preventing any 
reliable analysis. Alternatively, one can calculate 
k by fixing the mass of an hadron containing a 
quark of a given flavour to its experimental value 
using eq. ([!]). This is the Spectroscopy method. 
Axial Ward Identities (AWI) are obtained 
analogously by applying axial-vector transforma- 
tions. In this case, however, there is no partially 
conserved lattice axial current either for the Wil- 
son or SW-Clover action. Therefore, any lattice 
definition of the axial current get renormalized, 
even in the limit of vanishing bare quark masses. 
For degenerate quarks, close to the chiral limit 



(k = k c ) and neglecting terms of order O(a), one 
obtains the AWI § 
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where P = 927591 is the bare pseudoscalar den- 
sity, — <Z27p75<7i the bare axial current and 
Zp and Za are their renormalization constants. 
From (|4|), one can define the renormalized quark 
mass as 



TURd 



ZApa = Z A a (0|9 A |7r(p = 0)) 
Z P 2 Z P 2 (0\P\ir{p = 0)) 



(5) 



where \ir(p = 0)) is a pseudoscalar meson state 
containing a quark of hopping parameter k. For 
the SW-Clover improved case, the only differ- 
ence is that one has to consider the ratio of im- 
proved operators and take the symmetric deriva- 
tive. This is the Axial Ward Identity method 
which does not use xPT. 

From the VWI and the AWI some interesting 
properties of the renormalization constants can 
be obtained: 

1. Zy, where V M = 927/Wi, and Za are finite 
functions of g$ with Zy, Z A 7^ 1. 

2. Zp and Z$ are logarithmically divergent in 
/xa but Zs /Zp is a finite function of <?q only 
with Z s jZp + 1 for the Wilson and SW- 
Clover lattice actions. 

3. Evaluation of the renormalization con- 
stants. 

Matching between lattice and continuum is 
necessary to convert the lattice results to the con- 
tinuum renormalization scheme one has chosen 
to analyze the experimental data. The match- 
ing factors Z's are called lattice renormalization 
constants. At scales a -1 = 2 — 4 GcV of our 
simulations, one expects small non-perturbative 
effects on the Z's so that lattice (LPT) and con- 
tinuum (CPT) perturbation theory may be used 
to calculate them. In almost all cases, however, 
the Z's are known only to one-loop and, further, 
is well known that LPT corrections are large due 
to tadpole diagrams. It is claimed that the so- 
called Boosted perturbation theory (BPT) can re- 



3 



sum tadpole diagrams by replacing the bare lat- 
tice coupling constant with a non-perturbatively 
renormalized one, ap — cjq/(Att)P, where P is 
the plaquette expectation value, so that the ob- 
servables calculated with this coupling are closer 
to their non-perturbative counterparts. Unfor- 
tunately, data show that even after the imple- 
mentation of the BPT, 0(a 2 s ) corrections can be 
large, of the order of 10 - 50%, depending on the 
quantity at hand. Alternatively, one can use the 
non-perturbative renormalization technique pro- 
posed by Martinelli e£ aZ The idea is to im- 
pose renormalization conditions directly on quark 
and gluon Green functions in a fixed gauge with 
given off-shell external states with large virtu- 
alities /i 2 , mimicking the continuum renormal- 
ization procedure but performing all computa- 
tions non-perturbatively on the lattice. This 
procedure defines the Regularization Independent 
scheme (RI), also called MOM. For example, the 
renormalized scalar density in the RI scheme is 
P RI 0) = Z P (fj,,a)P(a) where Z P is fixed by 
imposing the renormalization condition 

Zp{^,a) (p\P(a)\p) | p2= _ M 2 = (p\P(a)\p) (6) 

where (p\P(a) \p) Q is the tree level matrix element. 
Since (^|) can be evaluated non-perturbatively 
through lattice simulations, Zp is extracted with- 
out LPT. Notice that Zp is affected by possible 
0(a) effects and depends on the external states 
and on the gauge. This dependence, however, will 
cancel in the final results at a given order when 
combined with the continuum MS-R1 matching 
calculation. For this method, which defines the 
same renormalized operators in all regularization 
schemes, to work the following window in [i must 
exist 



(7) 



where the first term is required to avoid large 
higher-order corrections in CPT while the sec- 
ond one is necessary to avoid 0(a) effects. The 
evolution of a bilinear operator renormalized in 
the RI scheme O^(p) {T = V,A,P or S) can be 
expressed in the form 
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where 1 are the solutions of the RGE which are 
known up to N 2 LO ||. Therefore, we can define 
a renormalization group invariant quantity 
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which, up to higher order terms in CPT, should 
be independent of /i, in the region in which 
CPT is valid /i > 2 GeV, of the renormalization 
scheme, of the external states and gauge invari- 
ant. In fig.l and 2, we show our preliminary re- 
sults for the Z's for the unquenched Wilson ac- 
tion (see ref. 0). We can clearly see a window 
in fi where the Z's are scale independent. The 
corresponding Z's for the Wilson and SW-Clover 
actions in the quenched approximation are simi- 
lar, but somewhat less stable, and can be found 
in ref . || . We have used them to perform the first 
measurement of quark masses without LPT. 

4. The continuum MS quark masses. 

Since the MS-Kl matching factor depends on 
typical scales of order [i m a -1 « 2 — 4 GeV, 
it can be calculated in continuum perturbation 
theory only. We consider two cases: 
Perturbative calculation of the Z's, i.e. one 
uses both BPT and CPT at NLO 
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where U^ s is the evolutor operator for the quark 
mass at NLO in th e co ntinuu m. T he re normaliza- 
tion constants \ jZf s and Zf s jZf s are known 
to 0(us) only (see ref. Q for details). 
Non Perturbative calculation of the Z's, 
i.e. one uses CPT but the Z's are computed non- 
perturbatively as explained in section 3. In this 
case, a N 2 LO matching can be performed 



m MS (ri = U% s W) 
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where and are M5*-RI matching con- 
stants (see ref.fl for details). 
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Figure 1. Zy and Za for the unquenched Wilson 
action as a function of fi at n sea = 0.1580. 
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Figure 2. Zp and Zg for the unquenched Wilson 
action as a function of /z at n sea = 0.1580. 



5. Results. 

We have analyzed about 1000 quenched con- 
figurations for the Wilson action and about 2000 
quenched configurations for the SW-Clover ac- 
tion in a variety of physical volumes (ranging 
from 16 3 x 32 to 24 3 x 64) and lattice spacings 
corresponding to f3 = 6.0, 6.2 and 6.4. We re- 
fer the reader to the ref. || for technical details. 
For the spectroscopy method, we find a reson- 
able agreement between the quark masses ob- 
tained with perturbative (PT) evaluated renor- 
malization constants and non-perturbative (NP) 
ones and are compatible with previous determi- 
nations. This not the case, however, for the AWI 
method where the PT results are lower than the 
NP ones by more than two standard deviations. 
Only when we use non-perturbative renormaliza- 
tion constants, the spectroscopy and AWI meth- 
ods give consistent values of the quark masses, as 
it should be. We think that the reason is that 
lattice perturbation theory fails to determine Zp 
even if the Boosted recipes are implemented. As 
for the charm quark mass, our results are less 
stables than for the light and strange quarks. We 
think that the reason is that there is a rather large 
0(ma) contamination in our data for the charm 
quark mass. Finally, we want to stress that we 
cannot extrapolate to the continuum limit, even 
though we have simulated at three values of [3, 



due to the small physical volume at (3 = 6.4. 
Since our results at (3 = 6.0 and (3 = 6.2 are defi- 
nitely stables within errors, i.e. we do not see any 
dependence on a in our results, we average the 
NP AWI results and obtain our best estimates: 

m^(2GeV) = 5.7 ±0.1 ±0.8 MeV 

mf^(2GeV) = 130 ± 2 ± 18 MeV 

mf^(2GeV) = 1662 ±30 ± 230 MeV (10) 

where the first error is statistical and the second 
is an estimate of the systematic error evaluated 
from the spread in the values of the quark masses 
extracted from different mesons, lattices and chi- 
ral extrapolation methods. 
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